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ABSTRACT 

NGC 7793 - S26 is an extended source (350 pc x 185 pc) previously studied in the radio, 
optical and x-ray domains. It has been identified as a micro-quasar which has inflated a super 
bubble. We used Integral Field Spectra from the Wide Field Spectrograph on the ANU 2.3 
m telescope to analyse spectra between 3600-7000 A. This allowed us to derive fluxes and 
line ratios for selected nebular lines. Applying radiative shock model diagnostics, we estimate 
shock velocities, densities, radiative ages and pressures across the object. We show that S26 
is just entering its radiative phase, and that the northern and western regions are dominated 
by partially-radiative shocks due to a lower density ISM in these directions. We determine a 
velocity of expansion along the jet of 330 km s _I , and a velocity of expansion of the bubble in 
the minor axis direction of 1 32 km s~ 1 . We determine the age of the structure to be 4. 1 x 1 5 yr, 
and the jet energy flux to be (4 - 10) x 10 40 erg s -1 The jet appears to be collimated within 
~ 0.25 deg, and to undergo very little precession. If the relativistic f3 ~ 1/3, then some 4 M 
of relativistic matter has already been processed through the jet. We conclude that the central 
object in S26 is probably a Black Hole with a mass typical of the ultra-luminous X-ray source 
population which is currently consuming a fairly massive companion through Roche Lobe 
accretion. 
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1 INTRODUCTION 

The Sculptor Group galaxy NGC 7793 is a SA(s)d spiral galaxy 
located at RA (J2000) 23h57m49.83s,Dec. -32d35m27.7s. It has a 
distance estimated from the tip of its Red-Giant Branch of 3.6 Mpc 
(Jacobs et al. 2009) or 3.44 Mpc estimated from its Cepheid stars 
(Pietrzynski 2010). The foreground Galactic Extinction E(B-V) is 
estimated at 0.019 magnitude (Schlegel et al. 1 19981) . In this paper, 
we adopt a distance of 3.5 Mpc. 

The elongated Hn region within it, N7793-S26, was discov- 
ered by Blair & Long ( 1997) using interference filters isolating the 
Hey + [N n] lines and the [S n] lines in order to identify Supernova 
Remnants (SNRs). They described the object as a extended oval 
emission region with a high [S n] to Ha ratio and a size of 185 x 350 
pc (based on the distance of 3.5 Mpc rather than the 3.38 Mpc used 
in their paper). The major axis of this nebula is oriented southeast 
to northwest direction. They noted there was no evidence of an inte- 
rior star cluster, but they suggested that multiple supernovae (SNe) 
may be in the process of creating a super bubble, which would ac- 
count for its unusually large size. 



Read & Pietsch (1999) reported X-ray observations of P8 co- 
incident within 3" to N7793-S26, best fitted by a low-temperature 
thermal model. They could not rule out the possibly of time vari- 
ability as seen in their Fig. 4 ROSAT light curves. Very Large Array 
(VLA) radio observatio ns at 6 and 20 cm followed, and were re- 
ported in Pannuti et al. J2002h . They noted the object is composed 
of both a bright point and extended source, designated together as 
S26. Their discussion of the object is quite extensive, noting it to be 
the only evolved stellar candidate in N7793 detected in the X-ray, 
optical and radio domains. The object was seen in Chandra and spa- 
tially r esolved as a triple source for the first time by Pakull & Grise 
( 2008), and was further studied by Paku ll et al. feplOh and Soria 
et al. d201dh . Later, Pannuti et al. ( feOllh also detected the X-ray 
counterpart to S26 in their analysis of the Chandra data obtained in 
NGC7793. 



Pakull & Grise ( 2008) noted that a significant fraction of 
Ultra-Luminous X-ray sources appear to be embedded in extended 
(~ lOOpc) bubbles of shocked gas, and they identified N7793-S26 
as belonging to this class of objects. Pakull et al. ( 120 lOh further 
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suggested that N7793-S26 is powered by black hole with a pair of 
jets, similar to SS433; but twice as large and more powerful. They 
suggest S26 has a structure similar to a Fanaroff-Riley type II ac- 
tive galaxy with a X-ray and optical core, X-ray hot spots and radio 
lobes encased in an optical and X-ray cocoon. The core can be fit 
by a power law in the 0.3-10 keV energy band, and has a luminos- 
ity of Lo.3-10 ~ 7 x 10 36 ergs s -1 . There are two X-ray bright hot 
spots on either side of the central source, which appear to be ther- 
mal in nature with a combined luminosity of L0 3-10 * 1.8 x 10 37 
ergs s _I , Using Hen narrow-band observations, they find signifi- 
cant emission at 4686 A (He 2+ — » He 1+ ) spatially distributed in a 
similar manner to Ha. A narrow slit spectrum taken near the core 
suggested an expanding shell with a maximum expansion velocity 
of ~ 250 km s _1 . 

Soria et al. feOld) used the Australia Telescope Compact Ar- 
ray (AT) to resolve the radio lobe structure and to map the spectral 
index in the radio cocoon of S26. The radio structure is strongly 
aligned to the axis defined by the three X-ray hotspots While the 
radio spectral index in the NW and SE lobes is approximately - 
0.7 to -0.6, it is flatter (-0.4 to 0.0) across most of the cocoon 
and becomes inverted at the base of the jets near the core (0.0 to 
+0.4J] suggesting a self-absorbed spectrum due to either free-free 
or Compton self-absorption. The total thermal energy in the bubble 
is estimated to be ~ 10 53 erg. 

In this paper, we present the results of integral field spec- 
troscopy in the wavelength region 3600-7 200A using the Wide 
Field Spectrograph (WiFeS) (Dopita et al. 12001 12010|) mounted 
on the 2.3m telescope located at the Siding Spring Observatory. 
This allows the shock velocity, pre-shock density, ram pressure and 
radiative age to be derived across the face of the nebula, in turn 
permitting a much more accurate determination of the fundamental 
properties of the micro-quasar which is driving these shocks. 



2 OBSERVATIONS AND ANALYSIS 

Integral Field Spectroscopy (IFS) on N7793-S26 was performed on 
13 December 2009 at Siding Springs Observatory using the 2.3 m 
Advanced Technology Telescope and its Wide Field Spectrograph 
(WiFeS). WiFeS provides a 25" x 38" field with 0.5" spatial sam- 
pling along each of 25" x 1" slits. The output format matches the 
4096 x 4096 pixel CCD detectors in both the red-side and blue- 
side cameras, and the instrument sensitivity is optimized for the 
blue and red sides of the spectrum reaching a throughput (top of 
atmosphere to back of camera) of 35%. We used the RT560 beam 
splitter. The blue spectral range covers 3200-5900A, at a spectral 
resolution of R ~ 3000 (~ 100 km s~'). In the red, the range is 
5300-9800A, and the resolution is R ~ 7000 (« 45 km s _I ). 

Three 1000 second exposures of S26 were taken under photo- 
metric conditions at PA (measured east of north) 145 degrees. The 
estimated seeing was excellent, sa 1.3". Standard and telluric stars, 
flat field, bias and wavelength calibration data frames were also ob- 
tained on the same evening. 

The data was reduced using the WiFeS data reduction pipeline 
based on NOAO (National Optical Astronomy Observatory) IRAF 
software. This data reduction package was developed from the 
Gemini IRAF package (McGreg or et al.. l2003h . Use of the pipeline 



consists of four primary tasks: wipes to set environment parame- 
ters, wftable to convert single extension FITS file formats to Multi- 
Extension FITS ones and create file lists used by subsequent steps, 
wfcal to process calibration frames including bias, flat-field, arc 
and wire; and wfreduce to apply calibration files and create data 
cubes for analysis. We used the standard star HD 26169 for flux 
calibration and the bright B-type star HIP 8352 for any necessary 
telluric corrections. 

In Fig. [T] we show a false color RGB image of N7793-S26 
in its proper orientation to the sky. There is intense [Sn] emis- 
sion in the southwestern hotspot, associated with the termination 
shock of the jet from the central object as seen in X-rays (Pakull 
et al. [2010|). The elliptical shell of the bubble is, by contrast, much 
brighter in [Om] and the [Om]/H/3 ratio increases as the Balmer 
line flux decreases. The shell brightness decreases going north on 
the eastern side, and is incomplete on its western side. This, to- 
gether with the asymmetry in the distance of the NW X-ray hotspot 
compared to the SW hotspot from the central X-ray source strongly 
suggest a systematic gradient in the density of the ISM, highest to- 
wards the SW side, and lower towards the NE. This is reminiscent 
of the gradient seen across the system SS433 - W50 in the Galaxy 
dZealev, Dopita & Malir]|l980h . 

The other emission regions appearing in this figure (the green 
region close to 23:58:00.5 -32:33:18 (J2000) and the turquoise re- 
gions in the general vicinity of 23:57:59.7 -32:40:00 (J2000) are 
H 11 regions which are unassociated with S26. These are useful in 
constraining the chemical abundances in the shock modelling de- 
scribed in the next section. 

The spectrum of the most intense 2" x 2" region is shown 
in Fig. [2] A classic shock-excited spectrum is revealed with 
both strong [Oi] and [Sn] emission relative to Ho- and velocity- 
broadened line profiles in all of the emission lines. This spec- 
trum can be compared with tha t of the emission region in SS443 
dZealev. Do pita & Maid dl98d) : fig 2), in which the velocity dis- 
persion is much lower, but the forbidden lines are much stronger 
with respect to Ha. The difference can be ascribed to two factors, 
first a much lower characteristic shock velocity in SS433 and sec- 
ond, lower chemical abundances in the ISM in S26. 

Using QFitsViewd we were able to create a Voronoi Tessella- 
tion diagram based on a 30: 1 Ha signal to noise (S/N) ratio using 
the red data cube as shown in Fig. [3] This was later used as a tem- 
plate to extract one dimensional spectra from NGC 7793-S26. After 
extraction, the IRAF task splot allowed the determination of emis- 
sion line flux and the velocity width based on Gaussian fits to the 
Ha line profile. Photometric measurement errors were calculated 
based on the root-mean-square noise estimated for each spectra; 
varying between an average of 2 percent for the stronger /16563 
line and 23 percent for the weaker /16364 line. 

In Table Q] we present flux densities (in units of 10~ 15 ergs 
cm -2 s -1 ) for selected emission lines corresponding to the Voronoi 
extraction by number listed in Fig. [3] The lines are labeled at 
rest wavelength and include H/3 014861), [Om] 0^4959,5007), 
[Oi] 01.16300,6364), [Nn] 0M6550, 6585), Ha (A6563) and [Sn] 
0U3716,3731). Also included are log ratios [Om]15007 / HyS, 
[Nn]i6585 / Ha, [S n]/U6716, 6731 / Ha and Ha surface bright- 
ness (ergs s _I cirr 2 sr~'). 



Spectral index here is denned as a, for S v cc v" where 5 „ is flux density 
and v, frequency. 



2 Written by Thomas Ott and freely available at: 
www.mpe.mpg.de/~ott/dpuser/index.html 
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Figure 1. RGB ([S n], Ha, [O m] /15007) false color image of NGC7793-S26. The y-axis has been rebinned to a final resolution of 1 X 1", and the coordinates 
given are for Epoch J2000. The purpose of this image is to map these three emission regions for visual inspection. However, the strength of the [S n] and Ha 
lines in the SE jet region is apparent, as is the relative strength of the [O m] lines in the fainter parts of the shell. The regions used for analysis are shown in 
Figure[3] and quantitive values for each of these emission lines in these regions can be found in Table[T] 



3 SHOCK ANALYSIS OF S26 

Although the SE region of S26 shows a very normal shock-excited 
spectrum, the relatively strong [O m] emission in the northern part 
of the shell is harder to explain. There are two possible explana- 
tions - either that this region is ionised by X-rays emanating from 
the central source, or that these regions contains only partially- 
radiative shocks which enhance the relative strength of the [O in] 
lines by cutting off the emission that would arise from the recom- 
bination zone of the shock. Such shocks were first considered by 
Dopita ( 119830 in the context of the supernova remnants IC 443 
and RCW 86, and the general behaviour of the line ratios of these 
"truncated" shocks were given. Later, Raymond et al. fl988l) made 
a more complete study, and investigated the effects of density and 
shock velocity. 

The possibility that the nebula is predominantly X-ray ion- 
ized can be readily eliminated on energetic grounds. The integrated 
luminosity of the X-ray sources in S26 is no greater than about 



3 x 10 37 erg s 1 from observations taken with the Chandra satel- 
lite (Pakull, Soria & Motch l2010Q . From the same reference, the 
H/3 luminosity of S26 is given as ~ 1 x 10 38 erg s~'. For compari- 
son, the sum of the H/3 flux in all the regions shown on Figure[3]is 
5.5 x 10 37 erg s _1 . However, we are probably missing some flux, so 
we adopt the (Pakull, Soria & Motch uOlOO value here. Models for 
photoionised nebulae with an abundance of Z = 0.5Z Q , appropriate 
to this region of NCG 7793, emit about 50 times the flux observed 
at H/3 (Lopez- Sanchez et al. 1201 20 . Therefore, the total nebular lu- 
minosity has to be of order 5 x 10 39 erg s - over 100 times the 
observed X-ray luminosity. 

Apart from these energetic constraints, if X-rays were domi- 
nating the excitation, the strength of the [O i] lines should be much 
greater in the SE portion of the nebula, since X-ray excited nebulae 
are characterised by extended warm, partially ionised tails in their 
Hn regions. This is expected on the basis of theoretical models as 
described in Dopita & Sutherland ( I2OO30 . and a fine example of 
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Figure 2. One-dimensional spectra of brightest 2" X 2" region in N7793-S26. Note the large velocity dispersion in the lines. Note also that the N7793-S26 
emission lines are redshifted (~ 228 km s -1 ) consistent with the reported redshift of N7793 (227 km s ; Koribalski et al., 2004). 




Figure 3. A portion of the elected Voronoi regions used for sampling overlaid onto part of the WiFeS Ha image. The 0.5 X 1 .0 arc sec. pixels of the WiFeS 
instrument have "stretched" this image in the horizontal direction compared with Figure 1. The regions 5, 6, 10, 12 and 13 are a complex of Hn regions, 
unassociated with S26. Apart from region 1 3 these are mostly located off the edge of this image in the lower left region. These H n regions were used to 
calibrate the chemical abundances appropriate to the analysis of S26. 



this phenomenon was observed by Pakull & Mirioni in the 

source Ho II X-l. 

We can therefore conclude that the spectral characteristic of 
the northern part of the nebula may be better explained by partially 
radiative shocks. This would be consistent with the evidence de- 
scribed above arguing for a lower pre-shock density in this part of 
the nebula. We now proceed to investigate this idea in more detail. 

3.1 The Shock Models 

The Hn regions identified in Figure [3] enable us to deter- 
mine the appropriate chemical abundance set. We have used the 
shock/photoionisation code Mappings Ills, an updated version of 
the code originally described in Sutherland & Dopit a i 19931) to 
generate the models described in Lopez-Sanchez et al. d2012h . We 
then used the line strengths reported in Table Q] and the modified 



strong -line technique also described in the Lopez-Sanchez et al. 
d2012l) to constrain the abundance in this region of NGC 7793; 
Z = 0.5 + 0.2Z o . We therefore adopt Z = 0.5Z o for our subsequent 
analysis of S26. 

With this abundance set, we have run a grid of partially- 
radiative shock models. The shock velocity was varied from 100 - 
220 kirr 1 , and the shock age was constrained by the condition that 
log([Om]^5007/ HB) < 1.0. We have assumed that the magnetic 
field pressure and the gas pressure are in equilibrium in the pre- 
shock gas, with an assumed pre-shock temperature of 10000K. For 
a pre-shock hydrogen density of 1 .0 this corresponds to a transverse 
component of the magnetic field of B = 1.68/iG. 

For all models, we have assumed full pre-ionisation of the gas 
entering the radiative shock. This assumption will be only valid for 
shock velocities in excess of 100km s -1 , so this is why we have 
restricted the minimum shock velocity to 100km s . 
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The key parameters and line fluxes of the models is given in 
Table [2] Column (1) gives, for each value of the shock velocity the 
corresponding shock age, expressed in terms of log(nf), where n is 
the pre-shock hydrogen density, and t is the shock age. All shock 
models with a given log(nf) will have very similar emission spec- 
tra, modulo collisional de-excitation effects in the density-sensitive 
lines such as [Sn]]6717,3lA. 

Column (2) of Table [2] gives the mean ion- weighted electron 
density enhancement factor in the zone emitting [S n], This ratio of 
post-shock to pre-shock electron density increases to a maximum 
as the gas is compressed by radiative cooling in the post-shock re- 
gion, before decreasing again as the shocked gas recombines, and 
the electron density decreases. Note that the computed compres- 
sion factor is larger for younder shocks. This apparently anoma- 
lous result is because the [S n] emission comes mostly from the 
peak of the electron density in the partially-radiative shocks, while 
in the fully radiative shocks the mean [S n] region electron density 
is much lower due to the recombination zone in the magnetically 
supported tail of the shock. The maximum gas density compres- 
sion factor is higher, and can be approximated by p\/po = 2 1/2 Ma, 
where Ma is the Alven Mach Number of the shock. 

Column (3) gives the surface brightness of the shock in Hy6 
for a pre-shock hydrogen density of 1.0 cm 5 . This will scale in 
proportion to the pre-shock density. 

The remaining columns of Table [2] give the strengths of some 
important emission lines with respect to H/?. These clearly show 
that the effect of decreasing shock age is to strengthen [Om] 
weaken [O i] and [S n] while the strength of [N n] is only slightly 
affected. 

3.2 Comparison with Observations 

In Figure [4] we show the computed line ratios compared with the 
observed line ratios across the face of S26. The observations imply 
shock velocities in the computed range of the models, 100 - 220 
km s -1 . The finite age models correctly predict the observed cor- 
relation between the (stronger) [Oiii]/H/3 ratios and the (weaker) 
[Oi]/Ha and [Sn]/Ha ratios. As expected there is little correla- 
tion between the [O m]/H/3 and the [N n]/Ha ratios. Qualitatively 
at least, the interpretation of the data as indicating the presence of 
partially-radiative shocks is confirmed. 

The analysis can be made more quantitative by comparing the 
surface brightness in Ha with the observed [O in]/Hy3 ratio. Models 
predict that for a given shock velocity as the shock becomes less 
radiative, the [O in]/Hy6 ratio should increase while the Ha surface 
brightness decreases. On the other hand, an increasing shock ve- 
locity leads to a higher Ha surface brightness as the mechanical 
luminosity of the shock increases, but the models show little varia- 
tion in the [O ni]/Hy3 ratio with velocity. 

In principle, these facts can be used to independently estimate 
both the shock velocity, V, and the radiative age log(/if)- However, 
we must first remember that in S26, we are looking through both 
the front and back of the shell, so we see two shocks along the line 
of sight. With this correction to the predicted flux we obtain the 
diagnostics shown in Figure [5] Here, the measured FWHM of the 
observed Ha line has been corrected for the instrumental resolution 
(50km s -1 ) by subtracting it in quadrature. Figure [5] implies that 
the lowest shock velocities are of order 80km s -1 , while the fastest 
shocks have V s ~ 220km s -1 . 

The shock velocities derived by this technique can be com- 
pared with the measured FWHM of the Ha line at the correspond- 
ing spatial position. We should not expect a one-on-one correlation 




FWHM (H-Alpha) (km/s) 

Figure 6. The correlation between the measured FWHM of the Ha line and 
the shock velocity inferred from Figure[3] The straight line would represent 
a one to one correlation. The existence of a correlation in this figure helps 
confirm that the estimated shock velocities are more or less correct. 

between these two values, since the FWHM does not necessarily 
reflect the actual shock velocity. It will be maximised when we 
observe two shocks moving in opposite directions, and minimised 
when we observe a shock travelling transversely to the line of sight. 
Nonetheless, since the radiative shock in general are observed mov- 
ing into denser clouds, and given the fractal distribution of these we 
might expect some correlation between FWHM and inferred shock 
velocity. Indeed there is, as is shown in Figure |6] We can conclude 
that the partially-radiative shock model provides a good description 
for the velocity width observations of S26. 

3.3 Constraints from He II 

Pakull et al. d2010h and Soria et al. < l2O10h claim (based on narrow- 
band observations) that there is a strong, extended Hell 4686 neb- 
ula, of similar size as the Balmer and [[Om]] nebula, and their es- 
timated He II 4686/H/3 flux ratio is ~ 0.1. We do not have data 
of similar sensitivity in the blue part of the spectrum, but an inte- 
gration across the nebula gives a flux ratio of He II 4686/Hy3 flux 
= 0.11 ± 0.033, in good agreement with the imaging data quoted 
above. 

Such a value is consistent wit h the f ast shock velocity, V s ~ 
250km s _I , inferred by Pakull et al. | |2010|) cf. Dopita & Sutherland 
Il996al. fi gure 6. The question we need to resolve here is whether 
such He II 4686/H/3 flux ratios are also consistent with the partially 
radiative but slower shocks we infer. 

Certainly, strong Hell is pr oduced in the precursor z ones o f 
fast shocks (Dopita & Sufherland [T996al . 1 1 996bl, Allen et al. l2008h . 
However, it must be remembered that the extent of the precursors 
is much greater than the shock itself, so any He II emission will 
be more spatially extended and more featureless than the shocked 
shell. 

Surprisingly, the slow shocks are also quite effective in pro- 
ducing He II emission (see Table [2] col. 4). Fully-radiative shocks 
produce He II 4686/H/3 flux ratios in the range 0.048 - 0.094 with 
a peak at about V s ~ 120km s -1 . This is consistent with the limit to 
the He II 4686/H/? flux ratio that can be placed upon the spectrum 
shown in Figure 2, taking into account the S/N of this spectrum. 
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Figure 4. A comparison of the line ratios derived from the partially-radiative shock models described in the text with the observational line ratios in S26. 
For clarity, a number of models have been omitted. The Red line is an isochrone for \og{nf) = 11.3 the yellow line for log(«?) = 11.6 and the blue line for 
logta) = 12.6. These diagnostics are insensitive to age when log(o/) > 12.6. 



The Physical Parameters of the Micro-quasar S26 in the Sculptor Group Galaxy NGC 7793 7 




-5.5 -5.0 -4.5 -4.0 -3.5 

log F(Ha) (erg cm -2 s _1 ster 1 ) 



Figure 5. The predicted relationship between the surface brightness in Ha and the [O iii]/H/? ratio as a function of shock velocity and radiative age log(nf). 
compared with the observations. Taken at face value, this figure implies shock velocities in the range 80 - 220 kirT 1 and radiative ages log(«?) > 1 1 in S26. 



Much higher values can be found in partially-radiative shocks, 
with values of He II 4686/H/3 > 0. 1 being readily attainable. Since 
partially radiative shocks will generally be moving into regions of 
lower density, and observational test would be to see if the fainter 
regions of the shell are characterised by higher He II 4686/H/? flux 
ratios. Unfortunately the S/N in the blue region of the spectrum is 
lower than in the red, and a combined spectrum of the faint regions 
(66-86) is just insufficient to provide a definitive test, since it only 
constrains He II 4686/H/3 < 0.15. 

3.4 Derived Shock Parameters 

In most regions of S26, the measured [Sn] ratio AA613 1/6717 is 
clearly above the low density limit of 0.66! 3, even when the mea- 
surement errors are taken into account. From Figure [5] we have al- 
ready estimated the shock velocity, V s and the radiative age param- 
eter, log(rcf). Thus using the observed [Sn] /1/16731/6717 ratio to 
infer the [S 11] density, and the appropriate [S 11] compression factor 
from Table (2) we can infer the pre-shock hydrogen density in each 
region of the nebula. 

The derived shock parameters are listed in Tablef3] This gives: 

• Col (1): the region number identified in FigurefJ] 

• Col (2): the inferred shock velocity from Figure [5] 

3 This value follows from the quantum mechanical sum rule for collision 
strengths between a singlet and a multiplet state; see equation (3.6) of Do- 
pita & Sutherland J2003I) . since at low densities the line flux ratio is simply 
decided by the ratio of the collision strengths from the ground state, cor- 
rected by the ratio of energies associated with the individual lines 



• Col (3): the measured FWHM of the Hff line. 

• Col (4): the age parameter log(rcf) inferred from Figure[5] 

• Col (5): the log of the inferred electron density, n e , from the 
[S 11] /1/1673 1 16111 ratio, where values of 1 .00 indicate that the ratio 
is measured at its low density limit. 

• Col (6): the inferred pre-shock hydrogen density, no- 

• Col (7): the age of the shock derived from columns (4) and 
(6), and finally, 

• Col (8): the ram pressure of the shock, P/k, derived from 
columns (2) and (6). 

In Figure UJ we show the correlation between the inferred age 
of the shock and the density and ram pressure in the shocks. Note 
that the analysis is insensitive to shock ages > 10 13 s. Broadly 
speaking, there is an inverse correlation between the pre-shock den- 
sity and the shock age, and that the pressure in the jet and counter- 
jet regions is larger than the average in the bubble. The inverse cor- 
relation between pre-shock density and shock age is to be expected, 
since the low density regions can only become radiative over rela- 
tively long timescales. 



4 GLOBAL PARAMETERS OF S26 

The theory of relativistic jet-driven bubbles was developed in the 
context of the double lobe radio sources in active galactic nuclei 
(AGN ) over many years by Scheue r Blandford & Rees 

1 19741) , Ra wlings & Saunders dl99ll). Kais er & Alexander d 19991) , 
Begelman |l996) and Bicknell et al. dl997h . Zealey, Dopita & Ma- 
lin J 198Ch developed the theory to apply to the case of the micro- 
quasar SS433. In summary, the relativistic jet is shocked when it 
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Figure 7. The correlation between inferred age of the shock and either den- 
sity (upper panel) or the ram pressure (lower panel). The red points on the 
right represent points in the SE jet and the two red points to the left of these 
diagrams represent shocks in the NW counter-jet. Note the high pressures 
in both jets, and the comparative (radiative) youth of the counter-jet shocks. 



starts to burrow into the ISM and the overpressure of the backflow 
carves out a cavity or cocoon around the jet. The back-flowing gas 
then inflates a larger bubble between the two jets, which behaves 
like a mass-loss bubble. The detailed physics is given in Dopita & 
Sutherland d2003h . Define the r-axis to be in the direction along the 
jets, and the z-axis to be the perpendicular direction. If the mean 
pressure in the cocoon is P, then towards the head of the cocoon (in 
the r direction) the mean pressure is higher, because this is defined 
by the small area,A, over which the termination shock is jittering. 
Here we can take the pressure as a factor f times the average lobe 
pressure where, typically f ~ 2 - 10. The velocity of advance of 
the jet, and the velocity of the wall shocks surrounding the cavity 
excavated by the jet are then given by: 



dr _ 
dt 

and 

dz = 
dt 



fit, 



pcA 



1/2 



1/2 



(1) 



(2) 



respectively where p is the density of the pre-shock ISM surround- 
ing the jet head, and Ej is the mechanical energy delivered to a 
single jet moving at a velocity /3c, and having an area A. 

For the larger bubble between the jets, Zealey, Dopita & Malin 

dl980h showed that the mass-loss bubble formula for the bubble ra- 
dius along the minor axis, z, given by Castor et al dl975T) or Weaver 
et al. dl977t) applies: 



,3/5 



(3) 



where Ej is now understood to be the combined energy flux of both 
jets, p ; is the density along the minor axis, and a is a constant. For 
a bubble containing thermal plasma with a filling factor 0, a = 
(125/154^-0) so that a l/5 = O.760~ 1/5 , while for a bubble entirely 
filled by relativistic plasma with y = 4/3, a = (125/2247T0) so that 

1.0) is more likely 



„l/5 



O.7O80~ 1/5 . The latter case (with , 



to be valid in the case of S26, since the extent of the non-thermal 
radio emission is quite similar to that of the optical shell^The age 
of the bubble, t, can then be directly inferred from the radius of 
the bubble, z, and the expansion velocity in the z-direction; r = 

(3/5)(z/v s ). 

The density, pressure and shock velocity in the cocoon around 
the SE jet can be obtained by averaging the values estimated for 
Spectrum numbers 24, 25, 28, 29, 32, 33, and 34. This gives 
v,(z) = 165km s -1 , P = 6.0 x l(r 9 dynes cnr 2 andp = 2.2 x 10~ 23 g 
cm -3 . The likely value of £ can be estimated by comparing the rate 
of advance of the jets given by equation [2] with the observations. 
With <T ~ 4 the rate of advance of the SE jet of 330 km s _I is repro- 
duced. Applying a similar analysis to the NW counter-jet (spectrum 
numbers 79 and 86) gives v s (z) = 130km s _1 , P = 2.5 x 10 -9 dynes 
cnr 2 , p = 1.5 x 10~ 23 g cnT 3 and a current rate of advance of the 
NW jet v,(r), of 260 km s _1 , which is also in good agreement to 
the observations. We may conclude that £ ~ 4. Finally, given that 
we would expect the the pressure at the termination shocks of the 
jets to fall off as the inverse square of the distance, the ratio of the 
measured pressures in the two jets should be in proportion to the 
inverse square of the distance. For S26 this implies a pressure ratio 
of ~ 2.2 which compares well to the observed ratio; 2.4. 

The velocity of expansion of the overall mass-loss bubble 
can be also estimated using the mean of spectra 66, 67, 70 and 
72. We find v s (z) = 132km s~\ P = 2.0 x Kr'dynes cnr 2 , 
p ; = 1.15 x 10~ 23 g cirT 3 . The ratio of v s (r)/v s (z) = 2.5 compares 
well with the ratio of the physical size in the r and z directions, 
R/Z =175pc / 92pc = 1.9, showing that the expansion velocity es- 
timates are reasonably consistent with each other. 

The age of S26 can now be estimated as t = (3/5)(z/v-) = 
1.3 + 0.15 x 10' 3 s or 4.1 ± 0. 5 x 10 5 yr. This is somewhat old er than 
the ~ 2 x 10 5 yr estimated bv lPakull. Soria & Motchl j20ld) thanks 
to the lower shock velocity estimated here. 

The jet energy flux can now be obtained from equation [3] We 



find E 



5.6^_x 10 4 ' 



identical to the value estimated 



e re 

by IPakull, Soria & Motchl d2010t) from the H/ 3 luminosity How - 
ever, they used the theory of fast shocks from I Allen et al.l (2008), 
whereas slow shock theory is more appropriate for the vast ma- 
jority of the shocks in S26. For these, the conversion factor be- 
tw een H0 (or Ha) fl ux and total luminosity is lower, and is given 
bv lRichet ail d2010l) . For shock velocities above 120 and less than 
200km s~', the conversion factor is virtually constant at about 80 



4 We are grateful to the referee for reminding us of the subtle distinction 
between the thermally and the relativistically filled cases. 
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(Ha to total luminosity). We therefore estimate the corrected total 
luminosity of the radiative shocks in S26 to be L rai = 2.6 x 10 40 
erg s _1 . Using the theory appropriate for a fully filled relativistic 
bubble in the fully radiative phase, Ej= (112/25)L rad = 1.0 x 10 41 
erg s~' (7.0 x 10 40 erg s" 1 if filled with thermal plasma). Combining 
these two approaches, we conclude that the jet energy flux lies in 
the range (4 — 10) x 10 40 erg s~', in excell ent agreement with that 
estimated bv lPakull. Soria & Motchl feOld) . 

The jet energy flux corresponds to an Eddington Luminosity 
for a 200 - 700M o Black Hole. However, it is comparable with 
the X-ray luminosities of the most hyper-luminous stellar sources 
(Swartz et al., 1201 ll) . strongly suggesting that S26 belongs to the 
same population of stellar Black Holes as the ULX sources. The 
difference is presumably that S26 is currently attempting to swal- 
low a stellar companion through Roche Lobe accretion so that the 
X-ray luminosity has been transformed to jet energy flux. 

By substituting the derived jet energy flux into equationQ]we 
can estimate the area A of the shocks. We do not have a good esti- 
mate of the relativistic f}, so retain this factor explicitly to estimate 
A ~ 0.38/T'pc 2 . If we assume /? to be similar to SS433 (J3 ~ 1/3), 
then the jet covers about l.lpc 2 , implying a collimation angle of 
about 0.25 degrees. 

If the X-ray hotspots are indeed marginally resolved, as 
claimed bv lPakull. Soria & Motchl (2010), then they cover ~ 20pc 
and subtend an angle of 6 degrees at the central source. Therefore, 
either the precession angle or "jittering" of the jets is within a cone 
opening angle of 6 degrees, or else the jet is de-collimated by the 
same amount. This latter condition would require a very low jet 
speed, fj ~ 0.02, of a jet velocity of ~ 7000kms~ 1 . Whilst not ex- 
cluded this seems unlikely, as it would require that a great deal 
of matter has been processed by the jet. This mass is given by 
M = r Ej l(J3cf. For /? = 1/3 we obtain M ~ 4M G , while for 
y8 = 0.02, the mass needed to be processed rises to over 1000M o . 
Clearly the jet energy flux cannot be maintained over the lifetime 
of the bubble unless the jets are fully relativistic. This also suggests 
that here we are dealing with a bubble filled with relativistic rather 
than thermal gas. The total energy injected to the bubble over the 
lifetime of the bubble is given by E wt =Ej t ~ 8 x 10 53 ergs. 



5 CONCLUSIONS 

We have shown that the radiative properties of the shell of NGC 
7793 - S26 can be understood in the context of a model in which 
the shell is still becoming fully-radiative in its NW portions, due to 
a large-scale density gradient in the ISM. Although the age of the 
oldest shocks cannot be reliably measured we can set an age limit 
comparable with the inferred age of the shell, 4.1 x 10 5 yr, and we 
can infer a jet advanc e velocity in the SE jet of 33 km s~ 1 . 

We concur with lPakull. Soria & Motor] d20ld) that S26 repre- 
sents the "missing link" between the ultra-luminous X-ray sources 
and their energetic bubbles, and the less luminous Galactic source 
SS433 and its nebula, W50. Indeed, the lower jet luminosity we in- 
fer here, Ej= (4- 10)x 10 40 erg s is comparable with t he X-ray lu- 
minos ities of the most hyper-luminous stellar sources ISwartz et aU 
J20T ]]). Unlike SS433, the jets in S26 undergo very little preces- 
sion; over not more than ±3 degrees. 

In conclusion, the central source of S26 is almost certainly a 
Black Hole of the typical mass which characterizes the ULX pop- 
ulation in a close binary with a more normal star. Given that S26 
has processed at least 4M G through the jets already, the companion 



is likely to be an intermediate-mass star which is in the process of 
being consumed by its unfriendly companion. 
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Table 1. Selected line fluxes (at rest wavelengths) and log ratios for N7793-S26 and the nearby Hn region to its south as shown in Fig. 3. Note in the case 
of [Sn]/ Ha, a ratio greater than 0.4 (the classical definition for a region dominated by shock emission; Lasker 1977; Mathewson et al. 1983, 1984, 1985) 
corresponds to a log value greater than -0.4. 
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The Physical Parameters of the Micro-quasar S26 in the Sculptor Group Galaxy NGC 7793 

Table 2. The finite age (partially radiative) shock models. All have a metallicity 0.5Z o . Line fluxes are given with respect to H/3 = 1.0. 



(1) (2) (3) (4) (5) (6) (7) (8) (9) 

log(«f) [Sn]Comp. log F Hp Hell [Om] [Oi] Hir [Nn] [Su] 

em's Factor erg craV'sr" 1 4686 5007 6300 6563 6584 6717+31 



V, = 100 kins" 



12.6021 


12.0 


-5.2206 


0.0481 


1.2090 


0.3544 


3.0970 


0.3728 


2.3861 


12.3010 


12.8 


-5.2251 


0.0486 


1.2210 


0.3581 


3.0860 


0.3767 


2.4111 


12.0000 


14.7 


-5.2410 


0.0504 


1.2670 


0.3707 


3.0770 


0.3908 


2.4910 


11.7782 


16.9 


-5.2640 


0.0531 


1.3360 


0.3767 


3.0760 


0.4119 


2.5240 


11.6021 


19.3 


-5.2991 


0.0574 


1.4480 


0.3606 


3.0820 


0.4452 


2.4720 


11.3010 


23.8 


-5.4150 


0.0746 


1.8910 


0.2869 


3.1210 


0.5510 


2.2865 


11.0000 


28.2 


-5.6755 


0.1335 


3.4450 


0.2082 


3.2780 


0.7975 


2.2044 



V s = 120 kms" 1 



12.6021 


15.8 


-5.0593 


0.0937 


1.4750 


0.3956 


3.0670 


0.4123 


2.5420 


12.3010 


17.0 


-5.0625 


0.0944 


1.4810 


0.3983 


3.0590 


0.4148 


2.5590 


12.0000 


18.8 


-5.0748 


0.0971 


1.5240 


0.4092 


3.0500 


0.4268 


2.6250 


11.7782 


21.9 


-5.0938 


0.1013 


1.5920 


0.4133 


3.0480 


0.4457 


2.6430 


11.6021 


24.3 


-5.1219 


0.1079 


1.6980 


0.3934 


3.0520 


0.4741 


2.5690 


11.3010 


30.6 


-5.2461 


0.1428 


2.2610 


0.2915 


3.0870 


0.5829 


2.2906 


11.0000 


35.5 


-5.4814 


0.2430 


3.8860 


0.2101 


3.2050 


0.7704 


2.1584 



V, = 140 kms" 1 



12.6021 


22.8 


-4.8116 


0.0656 


1.3780 


0.3733 


3.0160 


0.4738 


2.5500 


12.3010 


23.6 


-4.8132 


0.0659 


1.3830 


0.3747 


3.0110 


0.4755 


2.5600 


12.0000 


26.2 


-4.8215 


0.0671 


1.4100 


0.3815 


3.0030 


0.4847 


2.6050 


11.7782 


29.8 


-4.8375 


0.0695 


1.4630 


0.3810 


3.0010 


0.5026 


2.6060 


11.6021 


34.1 


-4.8710 


0.0749 


1.5800 


0.3463 


3.0040 


0.5391 


2.4950 


11.3010 


41.5 


-5.0405 


0.1099 


2.3340 


0.2346 


3.0390 


0.6505 


2.2247 


11.0000 


46.8 


-5.5030 


0.3107 


6.7670 


0.1524 


3.2840 


0.9212 


2.0544 



V s = 170 kms" 1 



12.6021 


30.3 


-4.6190 


0.0553 


1.0550 


0.3997 


2.9950 


0.5285 


2.7290 


12.3010 


31.4 


-4.6193 


0.0554 


1.0590 


0.4010 


2.9900 


0.5302 


2.7380 


12.0000 


35.1 


-4.6286 


0.0564 


1.0770 


0.4066 


2.9850 


0.5394 


2.7720 


11.7782 


40.0 


-4.6414 


0.0581 


1.1140 


0.3935 


2.9830 


0.5575 


2.7180 


11.6021 


46.7 


-4.6853 


0.0637 


1.2330 


0.3312 


2.9870 


0.6044 


2.5350 


11.3010 


56.8 


-5.0327 


0.1380 


2.7410 


0.1890 


3.0660 


0.7340 


2.1844 



V s = 200 kms" 
12.6021 
12.3010 
12.0000 
11.7782 



34.5 
37.2 
45.7 
62.1 



-4.4695 0.0589 

-4.4731 0.0594 

-4.4870 0.0608 

-4.5893 0.0728 



0.8596 
0.8650 
0.8948 
1.1320 



0.5061 
0.5102 
0.4868 
0.3173 



2.9940 
2.9880 
2.9860 
2.9970 



0.6038 2.8170 

0.6088 2.8410 

0.6281 2.7320 

0.7081 2.2227 



V, = 220 kms" 1 



12.6021 


37.8 


-4.3539 


0.0570 


0.7310 


0.5367 


2.9870 


0.6477 


2.6650 


12.4771 


38.9 


-4.3551 


0.0571 


0.7332 


0.5381 


2.9840 


0.6494 


2.6730 


12.3010 


42.6 


-4.3558 


0.0579 


0.7395 


0.5417 


2.9820 


0.6550 


2.6860 


12.1761 


47.0 


-4.3640 


0.0592 


0.7384 


0.5293 


2.9810 


0.5293 


2.6210 


12.0000 


59.1 


-4.3875 


0.0717 


0.7900 


0.4315 


2.9820 


0.6929 


2.2801 


11.9031 


72.9 


-4.5020 


0.1228 


1.0280 


0.2858 


2.9910 


0.7448 


1.8693 
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Table 3. The inferred shock parameters for each region of S26 defined by Figure[3] 



(1) (2) (3) (4) (5) (6) (7) (8) 

Spectrum V s FWHM log(«0 logn,, logn lo-gr log(i7*) 

Number Pons™') (kins" 1 ) (crrr 3 s) (cm™- 1 ) (cm™ 3 ) (yrs) (cnT 3 K) 



15 


80 


19 


80 


20 


130 


21 


80 


23 


80 


24 


140 


25 


140 


26 


80 


27 


1 10 


28 


170 


29 


170 


30 


120 


31 


100 


32 


170 


33 


170 


34 


140 


35 


130 


36 


170 


37 


180 


38 


130 


42 


120 


43 


190 


44 


200 


45 


140 


46 


120 


47 


170 


48 


200 


49 


170 


50 


170 


51 


140 


52 


140 


53 


170 


54 


200 


55 


200 


56 


110 


60 


1 10 


61 


140 


62 


170 


63 


100 


64 


140 


65 


120 


66 


140 


67 


80 


69 


80 


70 


140 


71 


120 


72 


140 


73 


110 


74 


100 


75 


110 


76 


130 


77 


100 


78 


100 


79 


120 


80 


100 


82 


80 


83 


80 


84 


120 


86 


80 



122 > 13.0 
97 > 13.0 

150 > 13.0 
129 12.6 
94 12.0 
147 > 13.0 
158 > 13.0 
97 12.0 

106 > 13.0 
153 12.6 

152 12.6 

119 > 13.0 
111 12.0 
166 12.6 

149 12.6 
111 12.6 
131 12.6 
165 11.7 

123 11.8 

97 11.5 
123 11.5 
173 11.8 

151 11.8 

98 11.4 
141 11.3 

184 11.6 

185 12.0 
151 11.6 
98 11.4 

192 11.5 
185 11.4 
157 11.6 
125 12.0 

123 12.0 
96 11.4 

193 11.6 
175 11.3 

120 11.4 
160 11.3 
155 11.3 
120 11.1 
157 11.2 

124 11.3 

107 11.3 

153 11.2 

125 11.1 
128 11.1 
116 11.1 
79 10.9 
110 11.0 
145 11.3 

108 11.1 
64 10.9 
151 11.2 
123 11.1 
68 11.2 
79 11.1 

150 11.2 
84 11.1 



2.20 


1.20 


2.30 


1.30 


2.76 


1.48 


2.95 


1.95 


2.00 


0.92 


2.40 


1.04 


2.60 


1.04 


2.40 


1.40 


2.4(1 


1.26 


2.40 


0.92 


2.4(1 


0.92 


2.70 


1.50 


2.00 


0.83 


2.40 


0.94 


2.40 


0.94 


2.05 


0.69 


2.00 


0.74 


1.80 


0.20 


2.20 


0.56 


1.00 


-0.52 


1 .60 


0.18 


2.20 


0.50 


1.80 


0.01 


2.20 


0.58 


2.23 


0.74 


2.00 


0.33 


2.00 


0.34 


2.00 


0.33 


2.35 


0.63 


1.80 


0.24 


2.30 


0.71 


1.80 


0.13 


2.00 


0.34 


2.00 


0.34 


2.35 


0.87 


2.00 


0.52 


2.35 


0.73 


2.00 


0.26 


1.10 


-0.27 


2.10 


0.48 


1.00 


-0.57 


1.95 


0.30 


2.30 


1.09 


2.60 


1.39 


2.45 


0.80 


2.15 


0.61 


2.40 


0.75 


2.45 


0.90 


1.00 


-0.46 


2.30 


0.78 


1.00 


41.54 


2.30 


0.87 


2.70 


1.23 


2.40 


0.92 


1.20 


-0.24 


1.00 


-0.10 


2.30 


1.20 


2.25 


0.75 


1.00 


-0.10 



4.30 


7.22 


4.20 


7.32 


4.02 


7.92 


3.15 


7.97 


3.58 


6.94 


4.46 


7.54 


4.46 


7.54 


3.10 


7.42 


4.24 


7.55 


4.18 


7.59 


4.18 


7.59 


4.00 


7.87 


3.67 


7.04 


4.16 


7.61 


4.16 


7.61 


4.41 


7.19 


4.36 


7.18 


4.00 


6.87 


3.74 


7.28 


4.52 


5.92 


3.82 


6.55 


3.80 


7.27 


4.29 


6.82 


3.32 


7.08 


3.06 


7.11 


3.77 


7.00 


4.16 


7.15 


3.77 


7.00 


3.27 


7.30 


3.71 


6.74 


3.19 


7.21 


3.97 


6.80 


4.16 


7.15 


4.16 


7.15 


3.03 


7.16 


3.58 


6.81 


3.07 


7.23 


3.64 


6.93 


4.07 


5.94 


3.27 


6.98 


4.17 


5.80 


3.35 


6.80 


2.71 


7.11 


2.36 


7.41 


2.90 


7.30 


2.99 


6.98 


2.85 


7.25 


2.70 


7.19 


3.86 


5.75 


2.72 


7.07 


4.29 


5.90 


2.73 


7.08 


2.17 


7.44 


2.78 


7.29 


3.84 


5.97 


3.80 


5.92 


2.40 


7.22 


2.95 


7.12 


3.70 


5.92 



